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Abstract 

In this paper, the fractional-order Burgers-Poisson equation is introduced by 
replacing the first-order time derivative by fractional derivative of order a. Both 
exact and approximate explicit solutions are obtained by employing homotopy 
perturbation method. The numerical results reveal that the proposed method 
is very effective and simple for handling fractional-order differential equations. 
Keywords: Fractional Burgers-Poisson equation; Homotopy perturbation 
method; Fractional derivative; Symbolic computation 

1. Introduction 

In 2004, the Burgers-Poisson (BP) equation has firstly been proposed to 
describe the unidirectional propagation of long waves in dispersive media [lj], 
denoted by 

U t + uu x = tp x , (1) 

(p xx = (p + u, (2) 

where (p and u depend on (t,x) 6 (0,oo) x i?, and subscripts denote partial 
derivatives. In order to well study BP equation, by applying 1 — 8% to ([]} and 
using ([2]) on the resulting right-hand side, we rewrite the BP system as a single 
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differential equation for u: 

ut - u xxt + u x + uu x = 3u x u xx + uu xxx . (3) 

Because BP equation is a shallow water equation modeling unidirectional water 
wave subject to weaker dispersive effects than the KdV equation. This means 
that BP equation is very important in the field of mathematical physics. In 
Ref. [l| , the authors turned out that BP equation features wave breaking in finite 
time, a local existence result for smooth solutions and a global existence result 
for weak entropy solutions were further proved. Later on, BP equation was also 



proposed by Fetecau and Levy |2j by Pade (2,2) approximation of the phase 
velocity that arises in the linear water wave theory. In Ref. [3], the authors 
used classical Lie method to construct group invariant solutions. Moreover, 
variational iteration method (VIM) was applied to study the numerical solutions 
of BP equation Ml. 

n 

Recently, fractional calculus has been extensively applied in many fields [5J. 
Many important phenomena are well described by fractional differential equa- 
tions in electromagnetics, acoustics, viscoelasticity, electrochemistry and mate- 
rial science. That is because of the fact that, a realistic modelling of a physical 
phenomenon having dependence not only at the time instant, but also the pre- 
vious time history can be successfully achieved by using fractional calculus. In 
particular, Wang 0, Q] employed the homotopy perturbation method (HPM) to 
solve the classical fractional KdV and KdV-Burgers equations, respectively. Mo- 
mani [§1 solved fractional KdV equation via the Adomian decomposition method 
(ADM) and Momani et al. [9] applied variational iteration method (VIM) to 
solve the space and time fractional KdV equation. In the above methods, HPM 
provides an effective procedure for explicit and numerical solutions of a wide and 
general class of differential systems representing real physical problems. Then 
HPM has been widely used by other authors [lfl as well as their referees to 
solve fractional order differential equations. 

In this paper, we introduce fractional-order into the BP equation ([3]) by 
replacing the first-order time derivative by fractional derivative of order a, then 



2 



we obtain the fractional Burgcrs-Poisson (fBP) equation of the form 

D"u - Dfu xx + u x + uu x - (3u x u xx + uu xxx ) = 0, i>0,0<a<l, (4) 

where a denotes the order of the fractional time-derivative in the Caputo sense. 
The function u(x, t) is assumed to be a causal function of time and space, i.e. 
vanishing for t < and x < 0. The general response expression contains a 
parameter describing the order of the fractional derivative that can be varied to 
obtain various responses. In the case of a = 1, the fractional Eq. Q reduces to 
the classical BP equation Furthermore, HPM will be employed to obtain 
both exact and approximate explicit solutions of fBP equation. 



2. Preliminaries 

We first give the definitions of fractional-order integration and fractional- 
order differentiation [12]. For the concept of fractional derivative, we will adopt 
Caputo's definition, which is a modification of the Riemann-Liouville definition 
and has the advantage of dealing properly with initial value problems. 

Definition 1. A real function f{t), t > 0, is said to be in the space C^, fi G R. 
if there exists a real number p > fi, such that f(t) = t p fi(t), where f\{t) (E 
C(0, oo), and it is said to be in the space C n if and only if f n £ C^, n G N. 

Definition 2. The Riemann-Liouville fractional integral operator of order a, 
J a , of a function f G C^, p, > — 1, is defined as 

J a f(t) = -L f (t- T) a -\f(T)dT, (a>0,t> 0). (5) 
r (") Jo 

Definition 3. The fractional derivative of f(t) in Caputo's sense is defined as 



1 f* / (m) (r) 

r(m - a) J Q (t - t) c 



DVit) - WZT—^ I v. , a \L m dr, (6) 



where m — 1 < a < m, m G N, t > 0, / G . 
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Definition 4. For m to be the smallest integer that exceeds a, the Caputo time- 
fractional derivative operator of order a > is defined as 

I — gt m i if a = m £ N, 

Lemma 1. If m — 1 < a < m, m e iV, and / € C™, /i > —1, i/ien 

m — 1 

J<*D"f(x)=f(x)-J2f(0 + )^- (8) 

fe=0 

3. Application of HPM to fBP equation 

Homotopy perturbation method (HPM) is first proposed by He [ijj]. It is 
a powerful mathematic tool to solve nonlinear problems, especially engineering 
problems. To further complement HPM, He has developed it Q, [ij]] recently. 

In what follow, the HPM is used to study the fBP equation @ with the 
initial condition 

u(x i 0)=x. (9) 

The exact solution of BP equation ([3]) , the special case of fBP equation (|4| when 
a = 1, is given by classical Lie method (see Ref.[^|) 

„(*,*) = I±£-i. (io) 

In view of HPM, the homotopy is constructed as following 

(1 - p)D^u + p(D?u - Dfu xx + u x + uu x - (3u x u xx + uu xxx )) = 0, (11) 

or 

D" u + p(-D"u xx + u x + uu x - (3u x ii xx + uu xxx )) = 0, (12) 

where p G [0, 1] is an embedding parameter. By utilizing the parameter p, the 
solution u{x, t) can expanded in the following form 

u(x,t) — u (x,t) +pui(x,t) + p 2 u 2 (x,t) + p 3 u 3 (x,t) + ■ ■ ■ . (13) 

Setting p — 1 gives the approximate solution 

u(x, t) = u {x 7 1) + ui(x, t) + u 2 (x, t) + us{x, t) H . (14) 
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Now, substituting (|14l) into (|T2"j) . and equating the terms with the identical 
powers of p, gives 

p° : DfuQ = 0, uq(x, 0) = a;; 

p 1 : Dfui - Dfu 0xx + u 0x + u u Qx - 3u Qx u 0xx - u u 0xxx = 0, u\(x, 0) = 0; 
p 2 : Dfu 2 - Dfm xx + ui x + (u ui x + uiu 0x ) - 3(u 0x u lxx + ui x u 0xx ) 

-(u Q ui xxx + uiu 0xxx ) = 0, u 2 {x, 0) = 0; 
p 3 : Dfu z - Dfu 2xx + u 2x + (u u 2x + uiu lx + u 2 u 0x ) 

—3(uq x u 2xx + u\ x u\ xx + u 2x uq xx ) — (uqu 2xxx + u\U\ xxx + u 2 uo xxx ) = 0, 

ua(x, 0) = 0; 

(15) 

Applying the operator J a on both sides of equations in ([15)) and utilizing Lemma 
[1] , yields 

uo(x,t) = x; 

ui(x,t) = — r(i+ Q )^ Q i 

„ (rr f\ - (l+^)(4r(l+a) 2 +r(2+q)) , 3a 
U3^, I)— T(T+a^ns+a) 



Since 

j'-i i-i 

Dt u j ~ + + X) u i u (j-i-l)x — 3 X Ui x U(j_i_x)xx 

i=0 i=0 nf\ 



j-1 

^i^(j — i—l)xxx 



z=0 

and Wj(a;,0) = 0, the other arbitrary Uj(x,t) (j > 4) can be calculated in 
the same manner by symbolic software programme Mathematica. If only the 
first four approximations of equation (|14[) are sufficient. Then the approximate 
explicit solution of equation (Q| will be expressed as 

( t \- 1 + x f a i 2 ( 1 + ^ 2a (l+x)(4r(l + a ) 2 +r(2 + a)) 3a 

[ ' j r(i + a) + T(2 + a) r(i + «) 2 r(3 + «) 

(18) 

In the other word, the exact solution for fBP equation when a — 1 can be 
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obtained by HPM. In face, if a = 1, equation (fTBJ) will read 

u (x,t) = x; 

ui(x,t) = -(l+x)t; ^ 
u 2 (x,t) = (l + x)t 2 ; 
u 3 (x,t) = -(1 + x)t 3 . 
According to (fT7|) . Uj satisfies 

Uj{x,t) = + x)t j . (20) 

Then the exact solution will be expressed as 

- 1 +x 

u(x,t) = lim 2_\uj(x,t) = — — - 1, (21) 

n->oo £ — ' 1 + I 

3=0 

n n 

which is just the same as for the classical Lie method [3| and the VIM [4 ] . 



Table 1: Comparison with the exact solution and the four-term approximation solution. 



(x,t) 


^exact 


UHPM 


\Uexact — UhPm\ 


(2,0.3) 


0.528462 


0.526 


0.002462 


(2,0.35) 


0.481481 


0.45925 


0.022231 


(2,0.4) 


0.428571 


0.392 


0.036571 


(2,0.45) 


0.37931 


0.32275 


0.05656 


(2,0.5) 


0.333333 


0.25 


0.083333 


(0.9,0.2) 


0.583333 


0.5804 


0.002933 


(1.2,0.2) 


0.833333 


0.8304 


0.002933 


(1.5,0.2) 


1.08333 


1.08 


0.00333 


(1.8,0.2) 


1.33333 


1.3296 


0.00373 


(2,0.2) 


1.5 


1.496 


0.004 



4. Numerical simulation 

In order to illustrate the approximate solution is efficiency and accuracy, we 
will give explicit value of the parameters x,t. Then calculate the two solutions 
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and make a comparison between them. Taking x — 2 for different value of t and 
t = 0.2 for different value of x, we calculate the numerical solutions of the two 
solutions which are given in (|10l) and (jT5J) . We list their numerical solutions in 
Table 1. 

From the numerical solutions in the Table 1, it can be seen that at the same 
time x, the value of the approximates solutions and the exact solutions are quite 
close. Also, when the value of t decreases the approximate solutions are more 
and more closed to the exact solutions. This shows the approximate solution 
is efficiency. It is also suggested that HPM is a powerful method for solving 
fractional differential equation with fully nonlinear dispersion terms. 

5. Conclusions 

In this paper, both exact and approximate explicit solutions of fractional- 
order Burgers- Poisson equation are obtained by employing homotopy perturba- 
tion method. The comparison of numerical solution and exact solution demon- 
strates that the proposed method is very effective and simple for solving solu- 
tions of fractional differential equations. It should be pointed out that detailed 
studies of fractional-order Burgers-Poisson equation are only beginning. The 
periodical waves, peakons, fractional Hamilton structure and other properties 
are still open. We hope that this work is a step in this direction and some 
traditional analytic method for nonlinear differential equations of integer order 
can be extended to fractional-order equations. 
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